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Abstract. The model of the Belousov-Zhabotinsky reaction, so-called 
hodgepodge machine, is discussed in detail and compared with the ex¬ 
perimentally determined system trajectory. We show that many of the 
features observed in the experiment may be found at different internal 
excitation growth rates, the g/maxState ratios. The g/maxState ratio 
determines the structure of the limit set and defines also details of the 
system state-space trajectory. While the limit set identical to the exper¬ 
iment and the model may be found, there is not a single g/maxState 
level which defines the trajectory identical to the experiment. We also 
propose that there is an inherent experimental time unit defined by the 
time extent of the bottleneck chemical reaction process which defines the 
g/maxState ratio and the spatial element of the process. 

Keywords: Chemical self-organization, multilevel cellular automata, spi¬ 
ral formation 


1 Introduction 

Until recently, the reaction-diffusion simulations based on PDEs, which are cen¬ 
tral models of chemical processes, have been used extensively to explain the 
self-organizing Belousov-Zhabotinsky (B-Z) reaction [I]. This kind of simulation 
expects an instantaneous chemical change. However, in reality, any chemical re¬ 
action is a combination of the physical breaking of chemical bonds, splitting, 
and diffusion of new chemical moieties over a time span - a defined elemen¬ 
tary timestep needed for the progression of the spatial-limited chemical process. 
Therefore, for the description of the B-Z reaction, we chose a kind of cellular 
automaton - a hodgepodge machine, see e.g. [2] - for the time-spatially discrete 
simulation. 

In brief, we consider four processes in the hodgepodge machine algorithm 
which describe the behavior of the chemical reaction at the quantum (electron) 
scale [§]: 
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1. Deexcitation phase transition - bond breakage upon reaching the maximal 
level, 

2. spreading the ’’infection” from neighboring cells, 

3. excitation process which is simulated by the acquisition of cell levels (energy 
transfer) from the average of neighboring cells and, 

4. growth of the excitation inside the cell. 

Items 1-2 are the most speculative processes, but we have at least the experi¬ 
mental evidence for them. Further, we assume that processes 3 and 4 occur on 
the border and inside of the structural element, respectively. In this aspect, the 
hodgepodge machine may be considered as the most elementary example with 
(a) only forward (excitation) processes and (b) a linear increase of the process 
inside the spatial element and on its border. The automaton depends upon 4 
parameters, maxState , k\, k, 2 , g , as well as on initial (an ignition phase) and 
boundary conditions (usually periodic), cf. [3j. 

In our previous paper [3j, we described the rule for the ignition phase of the 
hodgepodge machine which led to the course and final state highly similar to 
that observed in the B-Z reaction [4|-[6]: 

’’Let fci and be a weight of the cell in the ignition state from the interval 
(0, maxState) and that in the ignitial state maxState, respectively. Then, if 
1/Jfci + 0/^2 < 0.5 or 0/fci + 1 / k -2 < 0.5, a limit set of waves and spirals is 
observed. In all other cases we observe filaments.” 

In this article, we continue our research on the hodgepodge machine, studying 
the influences of the parameter maxState setting the maximal state value, and 
the internal excitation growth rate g , i.e. the number of levels increased in one 
simulation step. More precisely, we aim to understand how the ratio g/maxState 
influences the processes under the items 3-4. 

2 Methods 

We used the Wilensky NetLogo hodgepodge model 8], which was modified in two 
ways to improve its similarity to the B-Z experiment J3,9 . The first change (the 
generation of the random-exponential noise) was introduced to ensure that the 
simulation started with a small number of ignition points. The second modifica¬ 
tion (the rounding of the weighted occupied neighboring cells) was implemented 
in order to allow us to start from these few points, where the original rule (the 
calculation in integers) did not. The simulation was performed on a canvas of 
1000 x 1000 cells to avoid non-idealities caused by the periodic boundary con¬ 
ditions. Here, we present behavior of the simulation for k\ = k 2 = 3. 

For comparison to the simulation, the B-Z reaction was performed as previ¬ 
ously published [10 with the only modification in that we shook (1400 rpm) a 
big Petri dish (0 200 mm). 

In the discussion, we used Wuensche’s terminology 7|. The final states of the 
B-Z reaction and hodgepodge simulation are considered to be limit sets, and the 
courses of the experiment and its simulation are called trajectories through the 
state space. 
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3 Results and Discussion 

3.1 Influence of the maximal number of cell states maxState on 
the limit sets 

Multilevel cellular automata are often analyzed for at most 8 cell states, e.g. 

7 . Here, we tested the modified Wilensky hodgepodge machine algorithm by 
changing the maximal number of the achievable cell states maxStates (Fig. [I]). 

The limit set - a mixture of waves and spirals - highly similar to those ob¬ 
served in the B-Z experiment was achieved at the gj max State ratio of 28:200 
(compare Figs. [Tjz and ??c with Fig. [I] in 210 s). We observed branched spirals 
there (Fig. Bo- At the same ratio and higher maxState value (Figs. 01), the 
spirals were even more smooth. To obtain these results, we kept the ratio con¬ 
stant and simplified it to 7:50. The simulation at the g/maxState ratio of 3:20 
brought up the limit set, which was also similar to the chemical experiment. The 
g / maxState ratio of 1:7 (Fig. Bo led to the formation of square spirals reported 
earlier 17, Tl] . The next ratio 28:200 (close to 1:7) gave results most similar to 
the experiment. The g/maxState ratio of 1:8 (Fig. [TJd) already showed spirals 
and waves like the B-Z experiment. For all comonnly studied ” small-number” 
ratios, the simulation and reaction in the limit sets differ. At maxState < 20 
the spiral evolutions were already qualitatively similar. 

Among the low-number maxState levels, the simulation at the 1:7 ratio was 
extraordinary due to slow evolution to the limit set. Reaching the limit set took 
much higher number of simulation steps (37,000) than in other cases (e.g., 2,500 
steps at g/maxState = 1:8). Indeed, 8 levels, which correspond to the number of 
pixels in the Moore neighborhood, already tend to make octagons (read below). 
Just one number below, 7 levels, could change the system’s behavior. However, 
further research needs to be done to understand this phenomenon. 

3.2 Temporarily organized structures in the initial phases and the 
limit sets at different g/max State ratios 

The temporarily organized structures on the state-space trajectory of the simu¬ 
lation are those on one of the paths of the discrete dynamic network |7 . In Figs. 
[2] and [3] we show the shape dependence of the temporarily organized structures 
in the initial phases and limit sets on the g value at the constant maxState value 
(i.e., on the decreasing g/maxState ratio). 

Various ignition points, which are assumed to be Gardens of Eden, gave 
octagonal structures (Fig. [2]), whose interiors were typical for the given Gardens 
of Eden. At the high g/maxState ratio, after the passage of early waves when the 
interior octagon became regular, a centrally symmetrical structures appeared in 
the centre. The shapes of such initial central structures are dependent on the 
value of the g/maxState ratio. We do not describe here the mechanism how 
these structures arise. 

The early square waves gradually broadened with decreasing 
g/maxState ratio (Fig. [2]). The lower the g/maxState, the more diffuse the 
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Fig. 1 . The influence of the number of cell states maxState on the limit sets at similar 
g/maxState ratios of 1:7 (a), 3:20 (b), 28:200 (c), and 280:2000 (d). Total overviews 
of the limit sets’ structures. 


structures appeared. At g/maxState < 50:2000 we obtained a diffusive central 
object surrounded by the spreading wavefronts. With g/maxState of 10:2000, 
the fuzzy diffusive wavefronts were becoming more sparse. At g/maxState = 
1:2000 the second wave stopped evolving and the circular dark central structure 
appeared inside the dense diffusion. 

Thus, the g/maxState ratio can explain the thickness of the early waves and 
determine the number of states achieved in one time interval. This suggests that 
the wave thickness and density can, to some extent, explain other phenomena 
such as the built-in local asymmetry of the space [3], which ignites all additional 
phenomena and determines the state space trajectory of the process to its limit 
sets. 

The limit sets (Fig. [3]) first emerge from a fully spiral form at g/maxState 
= 1000:2000. As g/maxState decrease down to 280:2000, the spirals more and 
more resemble those found in a natural chemical reaction |3 . The waves were 
broadening and by g/maxState = 150:2000 mature spirals could no longer be 
observed. Later, the organized structures were broken into a diffusive mixture of 
darker and lighter stains. At the limit g/maxState —► 0, we are likely to obtained 
a uniform intensity. It shows the unsuitability of the reaction-diffusion model for 
the description of the B-Z reaction 1 . 
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Fig. 2. The initial states for selected values of g at the constant maxState = 2000. 
Each of the scale bars corresponds to 100 px. As in the previous figures, the black color 
corresponds to the state level 0, white color to the maximal state value (set by the 
constant maxState). 
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Fig. 3. The limit sets for selected g at maxState = 2000. From the upper left to lower 
right image, the values of g are 1000, 800, 500, 280, 150, 100, 20, 10, and 1, respectively. 
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3.3 Evolution of the state space trajectory 

The main goal of the hodgepodge machine simulation is to test its consistency 
with the course of the B-Z reaction (an experimental trajectory). The segments 
of the reaction are depicted in Fig. [IJ 



Fig. 4. Segments of the experimentally determined trajectory of the Belousov- 
Zhabotinsky reaction. Numbers correspond to the time interval (in seconds) from the 
initial mixing of the reactants. 


Despite the similarity of the limit set between the simulation at g/maxState 
ratio = 280:2000 (Figs. [l]: and ??c) and the chemical experiment, the trajectories 
towards these limit sets were a bit different. In the case of the simulation (Fig. 
[5]) , after the ignition, the structures evolved in a sequence of square waves with 
round centers. After that, the first spiral doublet arose in the vicinity of two 
central objects such that the waves collided in one point. The doublet soon 
became surrounded by an elliptical wave thicker than the early square wave. 
This elliptical wave system had an expanding ’’diffusive” character until it was 
compressed by another system of elliptical waves which slowly evolved at other 
points. 
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In case of the model, early circular waves analogous to those in the experi¬ 
ment, were achieved only at very low g/maxState ratios, when the waves were 
broad and the spirals were not formed (e.g., Fig. [2] for g = 20). This suggests 
that in the early phase of the experimental trajectory, there is a second process 
which brings the evolution of a different trajectory. Later, when this process is 
exhausted, the evolution of the system follows the path described by the hodge¬ 
podge machine model. 

4 Conclusions 

This paper examines the suitability of the hodgepodge machine simulation for 
the description of the Belousov-Zhabotinsky self-organizing reaction. The reason 
for doing so is that the hodgepodge machine represents a time-spatial discretiza¬ 
tion. The time-spatial discretization is the idealized process, where the time 
element determines (a) a set of events achieved at the time within a spatial 
element (a cell) the g constant - and (b) the rest of events on the spatial ele¬ 
ment’s (cell’s) borders - the maxState constant. We may consider analogies for 
such behaviour in known energy transfer processes, i.e. the difference between 
resonance energy transfer and energy transfer in systems of overlapping orbitals. 

The analysis of these constants showed that the g/maxState ratio determined 
the thickness of waves and the course of the state space trajectory irrespective 
to the total maxState value. The decreasing g/maxState ratio led to the loss of 
the wave structure, which eventually resulted in a fully diffusive picture. There 
also exists a lower limit of maxState for the appearance successful course of the 
simulation. The higher the maxState value, the smoother the spirals were. 
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- The octagonal waves in the centre collapses and is replaced by the fractal structure. 
Step 99 - In the places of the collision of two central structures, early spiral structures 
arise. Step 249 - Spiral structures are being surrounded by the elliptical wavefront. 
Step 749 The simulation canvas is being filled by new spiral waves. Step 1499 
Structures become almost regular by "compression” of waves of various origins. 








10 


Dalibor Stys et al. 


6. Zhabotinsky, A.M.: Periodic Liquid Phase Reactions. Proc. Ac. Sci. USSR, 157, 
392-395 (1964) 

7. Wuensche, A.: Exploring Discrete Dynamics. Luniver Press (2011) 

8. Wilensky, U.: NetLogo B-Z Reaction Model, available at 

http://ccl.northwestern.edu/netlogo/models/B-ZReaction Center for Connected 
Learning and Computer-Based Modeling, Northwestern Institute on Complex 
Systems, Northwestern University, Evanston, IL (2003) 

9. http://www.biowes.org 

10. Zhyrova, A., Stys, D.: Construction of the Phenomenological Model of Belousov- 
Zhabotinsky Reaction State Trajectory, Int. J. Comp. Math., 91(1), 4-13 (2014) 

11. Fisch, R., Gravner, J., Griffeath, D.: Threshold-Range Scaling of Excitable Cellular 
Automata. Stat. Comput. 1, 23-39 (1991) 


